Elastic transport through danghng-bond sihcon wires on H passivated Si (100) 
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We evaluate the electron transmission through a dangling-bond wire on Si(100)-H (2x1). Finite 
wires are modelled by decoupling semi-infinite Si electrodes from the dangling-bond wire with pas- 
sivating H atoms. The calculations are performed using density functional theory in a non-periodic 
geometry along the conduction direction. We also use Wannier functions to analyze our results 
and to build an effective tight-binding Hamiltonian that gives us enhanced insight in the electron 
scattering processes. We evaluate the transmission to the different solutions that are possible for 
the dangling-bond wires: Jahn- Teller distorted ones, as well as antiferromagnetic and ferromagnetic 
ones. The discretization of the electronic structure of the wires due to their finite size leads to 
interesting transmission properties that are fingerprints of the wire nature. 

PACS numbers: 73.63.Nm,73.20.-r,75.70.-i 



I. INTRODUCTION 

The development of information technology has been 
pursued at a tremendous pace. Larger capacity memories 
and faster processors are obtained from the miniaturiza- 
tion of electronic devices. Nevertheless, this technological 
explosion that started in the second half of the 20^^ cen- 
tury will reach a limit when facing the nanometer scale. 
Indeed, at this size, the bulk properties of semiconductors 
are not accessible any more. Thus, the operating prin- 
ciples will vanish along their miniaturisation process. 1^ 
In the mid-1970s, an alternative was proposed that sin- 
gle molecules could perform the same basic functions of 
electronics devices that are traditionaly silicon-based fab- 
ricated.^ ^ For this purpose, organic molecules are candi- 
dates with great potential given the control on molecular 
design rending possible by chemical synthesis. In particu- 
lar, one can conceive and experiment on single molecules 
that switch from one state to another under the applica- 
tion of some external stimulus^ or even amplify a signal.'^ 
Binary data can also be encoded in a single molecule to 
build up complex circuit. Then, logic operations, either 
basic (AND, NOT, OR) or more elaborated can be per- 
formed at the molecular level, giving rise to molecular 
logic gates .I^M^ 

One interesting problem is to design the atomic scale 
circuit where a number of these molecular logic gates will 
be (i) assessed and (ii) interconnected to obtain more 
complex circuits. Well designed atomic scale circuits can 
also perform logic functions by themselves without the 
need of molecular logic gates. It has been shown that 
such atomic scale structures can be constructed using 
the scanning tunneling microscope (STM). It can for ex- 
ample selectively remove H atoms on the surface of H- 
passivated semiconductors to construct dangling bonds 
( DBs) lines Those DBs are introducing electronic 
states in the surface band gap of the passivated semi- 
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FIG. 1: Atomic structure of (a) the Si(100)-(2xl)-H surface, 
(b) the infinite ideal wire drawn along the y direction. H 
atoms are depicted in cyan, while Si atoms are depicted in 
red (surface dimers), yellow (others) and blue when holding 
a dangling-bond. 



conductor surface. An important example of such con- 
struction is the DB wire produced by the selective re- 
moval of hydrogen atoms from an H-passivated Si (001) 
surface along the Si dimer row (see Fig. [l]). This DB wire 
has a single dangling bond per Si atom, intuitively offer- 
ing a ID metallic band within the gap. The transport 
electronic properties of such a wire has been previously 
inspected by Extended Hiickel calculations and explained 
using a tight-binding model. Unfortunately, the later 
configuration is not stable and a Peierls distortion takes 
place, involving a metal-insulator transition. The un- 
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stable DB wire, referred to as the ideal wire in the fol- 
lowing, can also relax in two magnetic forms resulting 
from the antiferromagnetic and ferromagnetic coupling 
of adjacent DBs, respectively. The description of these 
different states has been the subject of intense activity 
thanks to density functional theory (DFT) based calcu- 
lations f^^''^ In particular, when dealing with finite-size 
wires, the magnetic solutions appear to be more stable 
than the distorted surface structure, referred to as the 
non-magnetic (NM) wire.^^ Although these relaxed 
periodic structures break the metallicity of the DB wire 
by opening a gap, it does not dispose of the possibility 
to find a pseudo ballistic or a tunnel transport regime 
through a finite length DB wire. 

In the present work we inspect the transport properties 
of finite-size DB wires (ideal, NM, AFM and FM) using 
the non-equilibrium Green's Function (NEGF) approach 
combined with ab initio DFT. In order to complete the 
first model proposed by Kawai et a/.,^^ this study is pre- 
ceded by a thorough analysis of the infinite ideal wire by 
means of DFT-based calculations and Maximally Local- 
ized Wannier Functions (MLWEs).*^ 



II. COMPUTATIONAL DETAILS 



First-principles calculations are based on density func- 
tional theory (DFT) as implemented in Siesta .I^SEI cal- 
culations have been carried out with the GGA func- 
tional in the PBE form,^^ Troullier-Martins pseudopo- 
tentials,!^ and a basis set of finite-range numerical pseu- 
doatomic orbit als for the valence wave functions. Struc- 
tures have been relaxed using a double-^ polarized basis 
sets,^ while the conductance has been computed from 
first-principles, using a single-^ polarized basis set, by 
means of the TranSiesta method,!^ within the non- 
equilibrium Green's function (NEGF) formalism. In all 
cases, an energy cutoff of 200 Ry for real-space mesh size 
has been used. 

In order to reach a better understanding of the trans- 
port properties of the DB silicon wires, subsequent trans- 
formation of the Kohn-Sham orbitals to Maximally Lo- 
calized Wannier Functions (MLWFs)l^has been applied. 
This scheme allows one to obtain a localized orthogonal 
basis sets. As a consequence MLWFs offer an extremely 
convenient way to translate the problem in terms of an 
orthogonal tight-binding approach. Plus, an ab initio 
evaluation of the on-site energies and hopping integrals 
becomes available. The numerical calculations of MLWF 
have been run with the Wannier90 code,!^ used as a 
post-processing tool of Siesta. The interface between 
both codes has been developed earlier by R. Korytar et 



III. IDEAL WIRE AND H-JUNCTIONS 

The fully hydrogenated Si(lOO) surface presents a 
(2x1) reconstruction with dimer rows formed along the y 
{i.e. [110]) direction (see Fig. [l](a)). Starting from this 
fully passivated surface, one can construct a DB wire by 
removing H atoms along the dimer rows direction (see 
Fig. [l](b)). Before its relaxation, this DB wire presents 
a metallic character. Although this ideal DB wire relaxes 
in either a distorted wire (NM wire) or a magnetic phase 
(AFM and FM wires), the understanding of its electronic 
transmission properties remains a capital point as other 
cases will be regarded with respect to this reference. 
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FIG. 2: Transmission as a function of energy through the 
ideal wire (black solid line) and 1 (dashed red liner), 2 (blue 
dotted line), 3 (green dotted and dashed line) and 4 (purple 
double dotted and dashed line) H- tunneling junctions. 

Figure |2] shows the transmission of an infinite ideal 
wire. As expected from its band structure, the trans- 
mission exhibits clear steps featuring the existence of 
two channels between 0.08 eV and 0.26 eV above the 
Fermi energy (^f), with only one remaining channel for 
E — Ep taken between —0.53 and 0.08 eV. This result 
slightly differs from the one previously obtained by Kawai 
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from Extended Hiickel calculations, where the en- 
ergy range of the one-channel location is as large as the 
two-channel one (~ 0.4 eV), the total no zero transmis- 
sion in the surface gap being the same in both calcula- 
tions. 

The intuitive picture is that a transmission through 
such these non-relaxed DB wire arises from the direct 
through space coupling between DBs. However, the sub- 
strate is thought to play an important role in the trans- 
port,^^ and a convenient way to get rid of the direct 
coupling is to insert an H-passivated dimer in between 
two semi-infinite ideal DB wires (see Fig. [3|. This pe- 
culiar structure will be considered in section [V| to get 
a better understanding of the transport and scattering 
mechanisms in these systems. More passivated dimers 
can be added in between the two ideal DB wires to track 
the tunnel decay through such an atomic scale surface 
tunnel junction. 
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FIG. 3: Example of a 2H-j unction. H atoms are depicted in 
cyan. Si atoms are depicted in red (surface dimers), yellow 
(others) and blue when holding a dangling-bond. 

The transmission through different length H-j unctions 
is presented Figure [2] In the case of H-j unctions longer 
than one passivated dimer, a tunneling regime is observed 
in the middle of the gap with a transmission T{E) de- 
creasing exponentially with the length of the H-j unction. 
For a chosen energy, the transmission T{E) can be de- 
scribed as: 

where d is the distance between semi-infinite ideal DB 
wires in angstrom and 7 the inverse tunnel decay length 
(A"-*^). At the Fermi energy, the inverse decay length 
is found to be 7 = 0.29 ; a value comparable to 
the one 0.22 calculated with the Extended Hiickel 
semi-empirical approximation.^^ 7 can be monitored with 



0.38^ 
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FIG. 4: Variation of the tunnel decay length ^(E) (A~^) as a 
function of electron incident energy. Dashed lines correspond 
to the energies of the steps in the transmission of the ideal 
DB wire. 

respect to the energy (see Fig. |4|. It lies between 0.24 
and 0.37 A~^, and exhibits two local minima (maxima 
for the decay length Fig. [I] curve) corresponding each to 
the center of the one channel and two channels energy 
range of the ideal DB wire transmission spectrum (see 
Fig.[2|. 

The situation is qualitatively different for a single H- 
passivated dimer inserted between the semi-infinite wires. 
Indeed, this specific junction removes only one of the 



channels and does not affect the second channel (active 
between 0.08 and 0.26 eV). 

At this point, the ab initio DFT-treatment, although 
powerful, is too global to give a clear picture of the 
transport mechanism through surface DB wires and H- 
j unctions. In this regard, the parametrized tight-binding 
method constitutes an elegant way to model the physics 
of transport. Thus, one would like to take the best of 
both worlds and associate the accuracy of DFT-based cal- 
culations with the clarity of a simple tight-binding basis 
set. This is made possible by the use of Wannier func- 
tions, which allows one to (i) identify the main couplings, 
(ii) extract ab initio evaluations of on-site energies and 
hoping integrals, and (iii) use them in a tight-binding 
calculation. 



IV. WANNIER FUNCTIONS' ANALYSIS 

There is a large freedom in the choice of the Wannier 
functions used to describe the DB wire, in the sense that 
they are built to describe a certain energy interval of 
the band structure of the DB wire, and, also, the points 
in space around which they will be localized have to be 
specified as an input for the algorithm. In a first step, we 
choose to assign one MLWF at each bond and one func- 
tion per DB along the wire. The chosen energy interval 
was -12 < £; - < 3 eV. A subset of the MLWFs 
obtained for the ideal DB wire is depicted on Figure [5] 




FIG. 5: Plot of a dangling bond-like (left), and a subsur- 
face (right) Maximally Localized Wannier Functions. These 
MLWF are obtained by assigning one center per bond and 
allows one to recover the whole band structure of the system. 



With the corresponding Hamiltonian (written in this 
MLWF basis set), we can study the main electronic cou- 
plings in order to rationalize the electronic properties of 
the DB wire. One important feature is that the direct 
through space electronic coupling between nearest neigh- 
bour DBs is extremely low: —0.06 eV, as compared to the 
coupling between a DB and its subsurface functions that 
can reach —1.70 eV. Hence, one cannot think of trans- 
port as occurring only through direct hoping between 
DBs since the substrate dominates transport along the 
wire. 

From these functions, one can interpolate the band 
structure. The agreement with the ab initio DFT-based 
band structure is excellent. As for the transport proper- 
ties, we are mainly interested in the metallic band that 
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FIG. 6: Interpolated metallic band of the ideal wire along the direction of the DB wire. Color stands for the overlap of the 
eigenstate with MLWF corresponding to (a) the dangling-bonds, and (b) the dangling-bonds with 8 subsurface functions (see 
Fig. |5]). The Fermi level lies at —4.25 eV. The DB alone cannot account for the whole metallic band. Indeed, the band is well 
described only when taking into account 8 other functions. 



appears at the Fermi energy resulting from the presence 
of the DBs. Projections of the metallic band on different 
sets of MLWFs are depicted in Figure [6j Here, the metal- 
lic band is correctly described when including a large set 
of subsurface basis functions. This indicates that the 
price to pay for using bond-like orbitals is that a larger 
number of them is required in order to properly model 
the transport properties of non relaxed DBs chain, i.e. 
to describe the metallic band accurately. 

Alternatively, one can obtain a ID tight-binding de- 
scription of the system using only one basis function 
by considering the MLWF using only the metallic band. 
This leads to the calculation of an 'effective' orbital (see 
Fig. [t]) that is also delocalized over the substrate. 
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FIG. 7: Plot of the 'effective' MLWF. This MLWF is obtained 
by assigning only one center per DB and by describing only 
the metallic band. It can be seen as a weighed 'mix' between 
functions depicted in Figure [5] 



The metallic band can then be computed from the 
MLWF Hamiltonian, improving the calculations by in- 
cluding couplings to 2, 3 or more neighbors. We found 
that the coupling with the fifth neighbor had to be in- 
cluded in the calculation in order to obtain a quantitative 
agreement with the DFT-based result. In the next sec- 
tion, the model based on this description of the wire will 
provide an interpretation of the scattering of electrons on 
a single H-j unction. 



MODELING THE SCATTERING ON A 
SINGLE H-JUNCTION 



As presented earlier (see Fig. [2|, the single H-j unction 
has a transmission of ~ 1 in the energy window where it 
used to be T{E) = 2 without any H atom. Furthermore, 
it quickly drops to values close to zero for lower ener- 
gies. To model the presence of a saturated DB w hich 
acts as a scattering center, the scattering state^^^^ will 
be studied where the asymptotic states are the Bloch 
states of a periodic DB wire. We will compare the ab ini- 
tio transmission eigenchannels with the solutions found 
when modeling the DB wire as a one-dimensional chain 
with one orbital per DB site. 

We start by characterizing the available asymptotic 
(Bloch) states at each energy. The Figure [s] shows the 
band structure with the corresponding right and left- 
going states. They are labeled channel 1 and 2 respec- 
tively in the following. It can be seen that depending on 
the energy window, there are two or just one channel. 
When visualizing these states, the change of the phase 
for the effective orbital, at each site, is determined by 
the corresponding k value of the channel. 




R. Going Mode 1 

R. Going Mode 2 

L. Going Mode 1 

L. Going Mode 2 



k (a ) 

FIG. 8: Band structure of the ideal DB wire along the direc- 
tion of the wire for positive and negative values of k. Right 
and left-going modes are indicated. Above the blue (dashed- 
double dotted line) there are two modes (channels) that give 
a maximum transmission of two, and bellow it only one. a is 
the distance between two DBs. 
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We move on now to the scattering states. For ideal 
wires without any scattering center, the scattering states 
win correspond to the Bloch states at each energy (see 
Fig. [9| . Both channels 1 and 2 show the expected phase 
modulation corresponding to k vectors with values of 
ki ^ — 7r/(2a) and k2 ^ (where a is the distance be- 
tween two DBs). 
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FIG. 10: Scattering of an incoming wave into the available 
outgoing modes (here for a case where there are two modes, 
both, to the left and to the right of the scatterer), weighed by 
the corresponding transmission and reflection coefficients. 




FIG. 9: (a) and (b) display the channels 1 and 2 at E — 
Ef = 0.8 eV, respectively, (c) displays the scattering state at 
E — Ef = 0. Color code: real, positive in blue; real, negative 
in red; imaginary, positive in silver; imaginary, negative in 
orange. 

In the presence of an H atom (saturated DB), the in- 
coming states will be scattered as shown in Figure p^O) For 
energies values with just one channel and for a left incom- 
ing state, the resulting scattering state will be a super- 
position of the incoming wave (characterized by a certain 
kinc) and the reflected wave (with k = —kinc) weighed 
by the (complex) coefficient r; and, to the right of the 
scatterer, the transmitted wave will have k = kinc^ and 
will be weighed by the transmission coefficient t. As an 
example, in the case of the scattering state at E — Ep = 
(see Fig. [9] (c)), for a right-going state coming from the 



left, k 

7/ 



-7r/(2a) and T ~ 0.1. It can be seen that 



the transmitted wave corresponds to an outgoing wave 
with an associated k ~ — 7r/(2a). Since the transmission 
is low, with r ~ 1, the imaginary parts of the incoming 
and reflected waves cancel, leaving only the real parts of 
the amplitudes. 

The visualization process gets more complicated in 
the energy window where there are two channels, since 
now the scattering can couple the two channels (see 
Fig. 10). Although the transmission and reflection coef- 



problem and model the system as a one dimensional 
chain, as explained above, using the parameters extracted 
from the MLWF transformation. This strategy allows us 
to: (i) explicitly obtain the transmission and reflection 
coefficients and hence the total transmission coefficient 
T(E) (that depends also on the group velocities); (ii) 
obtain the transmission function in less than a second 
- instead of hours, the typical time it takes to obtain 
the T(E) from first-principles for the geometries consid- 
ered here. This latter aspect is particularly important 
if one wants to model more complex arrangements in- 
volving more than one saturated DB in extended tunnel 
H-j unctions. 

The saturated DB is modeled as a defect^ in the sense 
that, at the saturated site, there is no available orbital for 
the electrons to be transfered through, as depicted in Fig- 
ure 11 Different strategies can be used in order to solve 



the transport problem. Here we have used the Multiband 
Quantum Transmitting Boundary Method, proposed by 
Liang et al.^^ 
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FIG. 11: One dimensional model for the DB wires. Site 
represents a saturated DB, modeled as a defect, i.e. no trans- 
fer is possible through this site. With respect to the ideal DB 
wire, the Hamiltonian matrix element that is most affected is 
the one here represented by /3, followed by a (onsite), and /3^ 

When modeling the changes in the tight-binding ma- 
trix elements due to the presence of the H atom, one 
can expect that the most affected matrix elements are 
the electronic couplings between the two DBs separated 
by the saturated DB, in the sense of having their effec- 



ficients could be explicitly obtained from first-principles, 
we can now use the one-dimensional character of the 



tive through bond interaction reduced. Figure [TT] shows 
a scheme of the one-dimensional model chain with one 
saturated DB (sites 1 and -1 are separated by a satu- 
rated DB). The matrix elements indicated in Fig. [TT] 
are obtained from a MLWF treatment similar to the one 
applied on an ideal DB wire, and given in Table [ij The 
corresponding transmission curve is computed and shows 
an excellent agreement with the ab initio result. Also, for 
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TABLE I: MLWF Hamiltonian matrix parameters that were 
changed when going from the ideal to a wire with one satu- 
rated DB (together with the introduction of a defect). All the 
values are given in eV and correspond to the transfer integral 
matrix element shown in Figure pT] 



Chain type 


/3 


a 


ideal 


-0.053 


-4.333 0.190 


IH 


-0.041 


-4.298 0.182 



the energy window where there are two channels, the de- 
composition of the transmission function in terms of the 
contribution of the scattering of each channel into each 
other shows that they tend to be equally scattered. 

The effect of the passivation of one DB is then to create 
a scattering defect that effectively mixes the two incom- 
ing Bloch waves of different k values. This first inspection 
of the scattering properties leaves us with a physical pic- 
ture for the electronic transport in these systems. Bloch 
waves made up by the effective orbital, that are scattered 
by a saturated DB as a defect in a linear chain. Due the 
long range electronic interactions through the substrate, 
this one-dimensional chain has a somehow unusual band 
structure, such that there can be two channels for certain 
energies that correspond to two Bloch states associated 
to two different k values. This degeneracy gives rise to 
interesting scattering properties. 



VI. DANGLING-BOND WIRES 

As previously stated, an ideal finite DB wire is relax- 
ing by performing either a Jahn- Teller distortion (the NM 
wire) or a spin polarization with the DBs coupled ferro- 
or antiferromagnetically (FM and AFM wires). The 
structures and the energetics of these different atomic 
surface structures have been previously described. 1^^^^ 
In the case of short DB-wires, the AFM state is the 
most stable. But the calculated energy difference with 
other solutions remains too small to ignore 

them.™ 

Therefore, in this section, the transport properties of the 
three NM, AFM and NM DB wires configurations are 
discussed. 

The system is divided into left and a right leads and a 
central scattering region, in our case the different pos- 
sible wires (Ideal, NM, AFM and FM)."^^ The self- 
consistent density matrix is converged in the scattering 
region, using the open boundary conditions imposed by 
the leads through their self energies, using the now stan- 
dard Green's functions-based method. In order to avoid 
the sensitive problem of the description of the interface 
between realistic electrodes and the DB-wires, we choose 
to take advantage of the metallic behavior of the ideal 
wire (see section III) and use them as the source for 



electrons to be transferred through the scattering region. 
Thus, in our calculations, the role of right and left leads 
are performed by semi-infinite ideal wires. 

Notice that the central scattering region include (i) a 




FIG. 12: Setup used for transport calculations. Semi-infinite 
ideal wires (orange) are used as electrodes. The DB-wire is 
decoupled from the electrodes by H-passivated dimers (gray) . 



finite DB-wire with a length going from 2 to 5 DBs and 
(ii) two H-passivated dimers placed at each end of this 
DB wire (see Fig. 12 ). Those H-j unctions are used here to 
decouple the central DB wire from the leads. As a result, 
one can converge the different solutions corresponding to 
the different states of the DB-wire, while the leads remain 
as ideal DB wires. 
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FIG. 13: Band structure (left) and transmission T{E) (right) 
of the 5-DB ideal wire. Bands in blue (solid lines) correspond 
to states purely localized on the wire, bands in black (dashed 
lines) to those localized on the electrodes. The bands in red 
(dotted lines) result from a mixing between the states of the 
electrodes and those of the wire. The bands of the electrodes 
appear shifted due to the periodic boundary conditions used 
to calculate the band structure. 

As one can then expect, the effect of the two H- 
passivated dimers is to lead to a confinement of the DB- 
wire states, thus a particle in a box-like behavior. If in- 
stead of using the open boundary conditions, one uses pe- 
riodic conditions, these confined states will appear as es- 
sentially flat bands within the bulk gap as can be seen on 
Figure [13] on the example of a 5 DB ideal wire (blue solid 
lines). The corresponding transmission exhibits peaks 
and dips that can be understood by examining the con- 
fined states. 

The three first states of an ideal 5-DB wire give clear 
resonance peaks in the transmission spectrum. They are 
within the energy range of the one-channel energy range 
of the ideal DB wire lead, i.e. for (-0.53 eV < E - Ep < 
0.08 eV). On the other hand, the states located in the 
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FIG. 14: Same as Fig.fTslfor the 5-DB NM wire. 



two-channel energy range result from a mixing between 
the states of the DB wire and the ones from the leads. 
As a result, the transmission spectrum does not display 
resonance peaks but an overall enhancement in the high 
energy range of the two-channel band. 




FIG. 15: Plot of wavefunction displaying contributions on the 
wire and on the electrodes in the 5-DB NM wire (red dotted 
lines on Fig. 14) at F. 



The T(E) spectral behavior is similar for a NM wire 
(see Fig. 14). In the case a 5-DB NM wire, the DB wire's 
states give also rise to T(E) resonances. The larger trans- 
mission is obtained for a DB-wire state coupled with a 
state of the leads (see Fig. [l5|. In the two-channel en- 
ergy range, the T(E) spectrum is the one observed in the 
case of H-junctions (see Fig. |2|. Indeed, as opposed to 
the 5-DB ideal wire, the transmission is not perturbed by 
wire's states that are shifted and give resonance peaks at 
higher energies. When going from a 2 to a 5-DB NM 
wire, the overall shape of the transmission spectrum re- 
mains the same but with a larger number of resonance 
peaks corresponding to an increasing number of DB in 
the wire (see Fig. 16).^^ 

In the same manner, the transmission through a 5-DB 
AFM wire exhibits resonances corresponding to states lo- 
calized within the wire (see Fig. 17). Nevertheless, these 
states are confined below — £^f = — 0.3eV and above 
0.3 eV, leaving the transmission largely undisturbed. The 
situation is the same whatever the length of the wire (see 
Fig. 18). Let us stress that the difference between the 



transmission of both spins states in the case of odd finite 
length DB wires is due to the monoreferencial nature 
of the DFT approach that cannot deliver proper spin 
states. We observe an exponential decay of the trans- 
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FIG. 16: Transmission as a function of energy through NM 
wires. 2, 3, 4, and 5-DB wires are depicted in (a), (b), (c) 
and (d), respectively. The transmission of the ideal wire is 
depicted in black. 
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FIG. 17: Band structure (left) and transmission T(E^ (right) 
of the 5-DB AFM wire. The states localized on the wire 
are depicted in blue (solid and dashed lines for majority and 
minority spins, respectively). States of the electrodes are de- 
picted in black dashed lines. The bands of the electrodes 
appear shifted due to the periodic boundary conditions used 
to calculate the band structures. 



mission away from the resonances. Thus, one can de- 
fine and evaluate, as in the H-j unction case, an inverse 
decay length for the AFM wire. The calculated value 
lAFM = 0.89 A-^ dit E - Ef = 0.15 eV is over three 
times larger than the one of the H-junctions at the same 
energy. 

The case of the 5-DB FM wire is extremely similar to 
the previous one (see Fig.[l9|. The main difference lies in 
the energy distribution of the DB wire's states depend- 
ing on their spin. Indeed, as expected, the majority spin 
states occupy a lower range of energy, whereas the mi- 
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FIG. 18: Transmission through AFM wires as a function of 
energy, (a) 2-DB length, (b) 3 DB, (c) 4 DB and (d) 5 DB. 
Majority and minority spins are depicted in soUd and dashed 
hne, respectively. The transmission of the ideal wire is de- 
picted in black. 
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FIG. 20: Transmission through FM wires as a function of 
energy. The different pannels and lines correspond to the 
same cases of Fig. [18] 
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FIG. 19: Band structure (left) and transmission T{E) (right) 
of the 5-DB FM wire. The states localized purely on the 
wire are depicted in blue (solid and dashed line for majority 
and minority spins, respectively). States of the electrodes 
are depicted in black dashed lines. Bands in red dotted lines 
result from a mixing between electrodes' and wire's states for 
the minority spin. The bands of the electrodes appear shifted 
due to the periodic boundary conditions used to calculate the 
band structures. 



nority spin states are found near the two-channel energy 
range. As a consequence, the transmission spectra are 
extremely different for both spins, whatever the length 
of the finite central DB wire (see Fig. 20). The energy 



localization of the resonances is proper to the spin, fea- 
turing a spin filter behavior. 

This systematic inspection reveals a common feature in 
the transmission of finite size DB wires due to an overall 



general behavior of T{E). It is due to the original shape 
of the ideal DB wire, used here as electrode. This smooth 
transmission function is then perturbed by multiple reso- 
nances arising from states of the wire confined in between 
passivated dimers. The positions of these peaks and dips 
is specific to the nature of the wire as we have just seen. 

Figures [lOj [18] and |20] show the electron transmission 
coming from an ideal DB wire, transmitting into another 
ideal DB wire at a given electron energy E. As such, it 
is not easy to deduce the actual electron current in the 
studied systems. We have then computed the electron 
current using the Landauer-Biittiker formula*^ 



2e r 



T{E,V)[fR{E)-fUE)]dE. (1) 



Where, T{E^ V) is the transmission function for an elec- 
tron of energy E when the bias between the two DB 
electrodes is V, and fniE) {fL{E)) is the right- (left-) 
electrode Fermi occupation function. We further simplify 
the current / calculation using the zero-bias tansmissions 
of Fig. [16] [18] and [20] Figure [2l] shows the computed I-V 
curves for 5-DB wires in the ideal, NM, AFM and FM 
configuration. As we can see, the electronic currents of 
actual physical wires are very reduced as compared to 
the current of the ideally undistorted non-magnetic wire. 

Our recent study of total energy and stability of DB 
wires^^ shows that the AFM and NM solutions are ther- 
modynamically coexisting, however their very different 
I-V characteristics would permit to identify the type of 
obtained DB wire. 
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FIG. 21: Calculated current (nA) with respect to the bias (V) 
for 5-DB wires. Ideal, NM, AFM and FM wires are displayed 
top to bottom. 



VII. POLARONIC EFFECTS IN TRANSPORT 

The above study only considers elastic transport. How- 
ever, the large electron- vibration coupling leading the DB 
wires to Jahn- Teller distortions would modify transport 
from our above description. Indeed, a DFT-based tight- 
binding study suggested that transport in infinite DB 
wires would take place in the form of small-polaron dif- 
fusion .I^^Ell Hole polarons are extended over 3-DB sites, 
while electron polarons present a similar confinement.^ 
We have estimated the extension of the induced electron 
polaron, by including an extra charge in the NM 5-DB 
wire. The most noticeable effect is that the distortion re- 
verts sign: DB sites moving away from the surface when 
before they were moving against the surface and vicev- 
ersa. But the degree of distortion remains the same. 
More importantly, the electron is not localized to only 
three sites, but it extends over the full wire, giving rise to 
no confinement except for the finite size of the DB wire. 
Hence, the existence of small polarons may be due to the 
actual extension of the DB wire implying that transport 
will take place by polaron hopping only in long wires. 

The study of the electronic ground state of the neg- 
atively charged 5-DB wire further shows that there is 
little difference in the electronic states of the neutral 
and charged wires. Indeed, one can recover the other 



wire's electronic structure just by shifting the Fermi 
level. In this way, the new Highest Occupied Molecular 
Orbital (HOMO) is basically the former Lowest Unoc- 
cupied Molecular Orbital (LUMO) of the neutral 5-DB 
wire. Therefore, despite the need of including electron- 
vibration coupling in the quantitative description of DB 
wire transport,^ our present results are qualitative trans- 
port descriptions because the electonic current is largely 
tunneled via electronic states closely resembling the neu- 
tral DB-wire states. 



VIII. CONCLUSION 

In this study, we investigated the electronic transport 
properties of dangling-bond (DB) silicon wires on H pas- 
sivated Si(lOO). Thanks to DFT calculations and its 
analysis by maximally localized Wannier functions, we 
have been able to rationalize the transport properties of 
these DB wires. As an example, we have shown the effect 
of a Hydrogen impurity as a source of electron scatter- 
ing in the wire. This study shows that transport mainly 
proceeds via subsurface atoms because the direct DB in- 
teraction is negligible. A Hydrogen impurity is then ef- 
ficient in interrupting the subsurface transport because 
it decouples the subsurface states from the surrounding 
DB. In this way, a Hydrogen impurity mixes and scatters 
the channels of the otherwise unperturbed wire. 

Different finite-size dangling-bond wires have been 
studied. We have considered unperturbed, Jahn- Teller 
distorted, antiferromagnetic and ferromagnetic wires 
containing 2, 3, 4 and 5 DB's. Each wire displays typ- 
ical and unique trends in the transmission, allowing us 
to characterize the nature of the wire from its transport 
properties. The size of the studied wires are in the order 
of the small-polaron extension of DB wires, producing 
a shifting of the electronic structure but qualitative sim- 
ilar features. Hence, we expect polaronic effects mainly 
to be important in long wires. The transport properties 
revealed in this study will permit to characterize the DB 
wires experimentally accessible.^ 
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